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ABSTRACT 

We present theoretical results of a piecewise isothermal shock wind model devised 
for predicting the luminosity and surface brightness profile of diffuse X-ray emissions 
primarily from the inner shocked downstream wind zone of a planetary nebula (PN) 
surrounded by self-similar shocked dense shell and outer slow AGB wind envelope 
involving self-gravity and compare/fit our computational model results with available 
observations of a few grossly spherical X-ray emitting PNe. Matching shocked piece- 
wise isothermal self-similar void (ISSV) solutions with self-gravity of Lou & Zhai (LZ) 
for the outer zone and a stationary isothermal fast tenuous wind with a reverse shock 
for the inner zone across an expanding contact discontinuity, we can consistently con- 
struct dynamic evolution models of PNe with diffuse X-ray emissions. On the basis 
of such a chosen dynamic wind interaction model, both X-ray luminosity and radial 
X-ray brightness profile are determined by three key parameters, namely the so-called 
X-ray parameter X, two radii i?,.^ and Rc of the reverse shock and the contact dis- 
continuity. We find that morphologies of X-ray emissions would appear in the forms 
of either a central luminous sphere or a bright ring embedded within optically bright 
shells. In contrast to previous adiabatic models, the X-ray brightness peaks around 
the reverse shock, instead of the contact discontinuity surface just inside the outer 
shocked dense shell. Diffuse X-ray emissions of a few observed PNe appear to support 
this wind-wind dynamic interaction scenario. 
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1 INTRODUCTION 

Planetary nebulae (PNe) are produced at the late stage 
of stellar evolution, when a star transits from the asymp- 
totic giant branch (AGB) phase where a star blows slow 
AGB dense wind to the later phase showing a central com- 
pact star (i.e. hot white dwarf gradually exposed). As the 
later tenuous fast wind chases and collides with an earlier 
circumstellar slow massive AGB wind envelope, a forward 
shock emerges and runs outwards into the slow AGB wind 
zone, while a reverse shock travels radially inwards relative 
to the inner fast wind; a contact discontinuity in between 
separates the forward and reverse shocks. Such dynamic 
shock interaction of wind-wind collision rearranges the dis- 
tribution of AGB wind envelope and eventually result in a 
quasi-spherical PN (e.g ., Kwok et al. I. 1978; V olk fc Kwok . 
1 19851 ; IChevalier I . ll997i : lBian fc Lou 1 . 12005 ; ,Yu et al. 1 . 12006 ; 
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I Lou fc Zhail . l2009l , - referred to as LZ hearafter). Ignoring 
the gravity effect completely. Chevalier (1997) presented an 
isothermal dynamic model for PNe and constructed spheri- 
cally symmetric global hydrodynamic solutions to describe 
the expansion of outer shocked shell with an inner contact 
discontinuity o f fas t wind moving at a constant speed in 
the inner zone. ILZ I advanced a piecewise isothermal model 
framework including the gas self-gravity and showed var- 
ious plausible self-similar dynamic behaviours of outer PN 
envelopes. The dynamic collision interaction of the inner fast 
tenuous wind with the outer slow AGB wind envelope is well 
taken into account by inclu ding forward and reverse shocks 
and contact discontinuity in lLZ I . As an illustrating example 
of PN a pplic ation, isothermal self-similar void (ISSV) solu- 
tions of ILZ I appear sensible in fitting the available data of 
the PN NGC 7662. 

Diffuse X-ray emissions due to spectral line cooling of 
highly ionized metals (e.g. Kastner et al. 2000; Yu et al. 
2009) and probably partial thermal bremsstrahlung from 
hot electrons reveal important inner morphological features 
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of PNe, i.e. the region enclosed within the dense shell. A 
powerful reverse shock can heat up the inner tenuous fast 
wind plasma on the downstream side withi n an estimated 
temperature range of ^ 10^ — 10^ K (e.g . Volk fc Kwok 
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20061 ) and create a reverse shock heated 'hot bubble' largely 
confined by the outer dense shell at an expanding contact 
discontinuity. For this overall physical scenario, properties of 
diffuse X-ray emissions are mainly determined by the tem- 
perature, density, velocity and mass loss rate of the shocked 
inner wind. As a result, it offers unique insight for the shap- 
ing processes of PNe and their central hot bubbles. 

In the past decade, advanced instruments in space al- 
low us to detect difi'use X-ray emissions from PNe. Together 
with complementary observations for dense shells of PNe in 
optical and infrared bands, we are now able to test the wind- 
wind dynamic shock interaction scenario on more specific 
terms. X-ray emissions from three PNe were detected with 
ROSAT spacecraft (e.g. Stute & Sahai 2006), and from nine 
PNe with Chandra and XMM-Newton X-ray Observatories 
(e.g. Kastner et al. 2000, 2001, 2003; Chu et al. 2001; Guer- 
rero et al. 2002, 2005; Sahai et al. 2003; Montez et al. 2005; 
Gruendl et al. 2006). With sufficiently high spatial and en- 
ergy resolutions as well as detector sensitivity, Chandra and 
XMM-Newton spacecraft observations allow us to study the 
extended spatially resolved X-ray structure of PNe, which 
provide an unprecedented opportunity to examine the wind- 
wind dynamic interaction scenario. 

Akashi et al. (2006) qualitatively calculated the mor- 
phology of X-ray emissions from PNe based on the adiabatic 
self-simil ar colliding wind solut ions for shocked fast wind 
flows of I Chevalier fc Imamura I (J_983,) without self-gravity 
and showed that the X-ray emission appears in the form 
of a narrow ring inside the optically bright part of a PN. 
This was suggested to be the case of NGC 40 by Montez et 
al. (2005). Akashi et al. (2007) investigated X-ray emissions 
from PNe by numerical simulation of adiabatic winds with- 
out gravity. They found that in order to explain the observed 
PNe X-ray properties, a rapidly decreasing mass loss rate is 
required. The adiabatic colliding wind models of Chevalier 
& Imamura (1983) and Akashi et al. (2007) both indicate 
that X-ray emissions become brightest very near the dense 
shell. Therefore, PNe with central luminous X-ray emission 
profiles cannot be accounted for by their adiabatic collid- 
ing wind models. Chevalier & Imamura (1983) also included 
an isothermal dense shell in their model, in addition to the 
innermost and outermost stationary isothermal winds. 

Recent observations reveal that diffuse X-ray emissions 
from PNe do generally lie within the interior of the optically 
bright shells of relatively dense ionized gas, such as PN BD 
-F30°3639 (e.g. Kastner et al. 2000; Yu et al. 2009), NGC 
40 (e.g. Montez et al. 2005) and NGC 2392 (e.g. Guerrero 
et al. 2005). From PN BD -H30°3639, ROSAT spacecraft 
instrument discovered a Gaussian-shape diffuse X-ray emis- 
sion source (e.g. Leahy et al. 2000). With a high spatial 
resolution of 0".49, Chandra spacecraft well resolved this 
compact X-ray emission source and found a quasi-spherical 
central bright diffuse X-ray emission blob embedded inside 
the optical bright shell, with a radius ~ 2". The X-ray pho- 
ton flux from the centre is ~ 10 times more intense than that 
from ~ 2" (e.g. Kastner et al 2000). Moreover, Chandra also 
found the X-ray bright centre was slightly concentrated to- 



wards the east side. Different from PN BD -(-30°3639, PN 
NGC 40 shows a distinct annulus X-ray emission e.g. Montez 
et al. 2005). For those PNe of bipolar morphologies that do 
not clearly present shell-like structures, their diffuse X-ray 
emissions are also from within the optically bright shells, 
such as NGC 7027 (e.g. Kastner et al. 2001), NGC 7009 
(e.g. Guerrero et al. 2002) and NGC 6543 (e.g. Chu et al. 
2001), NGC 7026 (e.g. Gruendl et al. 2006) and Menzel 3 
(e.g. Kastner et al. 2003). These X-ray observations bear 
strong implications that wind-wind dynamic interactions or 
jet- wind interactions are most likely responsible for the inner 
diffuse X-ray emissions from PNe (e.g. Kastner 2007). 

This paper establishes a quantitative connection of self- 
similar piecewise isothermal dynamic evolution model of 
PNe with diffuse X-ray emissions, especially the rad ial dis- 
tribution of surface brightness. On the basis of |LZ ! global 
ISSV model for the self-similar evolution of PNe, we show 
that three key model parameters determine the morphology 
for diffuse X-ray emissions from PNe and speciflcally exam- 
ine observations of PN NGC 40. In our piecewise isothermal 
model framework, either ring-like or central bright X-ray 
morphologies are possible because X-ray emission peaks near 
the reverse shock instead of the contact discontinuity. Our 
results lend support to the wind-wind dynamic interaction 
model for PNe in general. We also compute the luminosity, 
surface brightness of X-ray emission from the PN 7662 and 
emphasize that this PN should be a highly desirable target 
source for Chandra further X-ray observations in the near 
future. 



2 A SPHERICAL WIND-WIND DYNAMIC 
INTERACTION MODEL FOR PLANETARY 
NEBULAE WITH SELF-SIMILAR SHELLS 

In this section, we first describe our global piecewise isother- 
mal dynamic evolution model for quasi-spherical PNe. At 
the late stage of stellar evolution, a star enters the AGB 
phase. During this phase, a star suffers from a rather high 
rate of mass loss ~ 10~^Mq yr~^ and thus creates a fairly 
dense AGB win d shell expanding at a slow radial speed of ' 



10 km s (e.g. 



1 19831 : IChevaiier 



Kwok et al. 1. Il978l : [ciievaher fc Imamura I, 



19971 . ILZ f ) . The dense gas shell keeps ex- 



panding until the stellar hydrogen envelope has been almost 
depleted and a compact hot white dwarf gradually exposed 
at the centre. Meanwhile, a hot tenuous fast wind somehow 
emerges from the central white dwarf as indicated by ob- 
servations. This tenuous fast wind blows at a high speed of 
> 10^ km s"^ with a mass loss rate of ~ 10"^° - 10"'' M© 
yr"'^ (e.g. Cerruti-Sola & Perinotto 1985). This fast inner 
wind will not diffuse through the outer AGB slow wind, but 
smash onto the outer dense AGB wind shell (e.g. Kwok et 
al. 1978) and then give rise to a forward shock propagat- 
ing outwards into the slow AGB wind and a reverse shock 
propagating inwa rds in the fast wind comoving framework 
of reference (e.g. IChevalier fc Imamura I (|l983l ): Akashi et 
al. 2006). A contact discontinuity naturally forms when the 
gas pressures of fast wind and of AGB wind reach a dynamic 
balance with both density and temperature being different 
on the two adjacent sides. These two shocks, especially the 
forward shock, rearranges the distribution of gas ffow and 
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finally results in a dense shell surroundi ng the centra l com - 
pact star blowing an inner fast wind re.g. lKwok| [2000l:lLzTl. 

A quasi-spherical PN may be roughly divided into four 
regions (see figure 1 of Akashi et al. 2006) . The fast wind oc- 
cupies the inner region from the central compact white dwarf 
to the reverse shock front at radius Rrs- When propagating 
inwards relative to the fast wind, the reverse shock heats up 
and slows down the inner wind and creates a downstream 
wind region from the reverse shock front to the expand- 
ing contact discontinuity at radius Re- These two regions 
form the inner zone and are characterized by the fast wind 
mass loss rate, wind velocities, densities, temperatures and 
reverse shock speed. The shocked dense shell region stays 
between the contact discontinuity and the forward shock at 
radius Rs- The fourth region, or the outer AGB slow wind 
envelope, expands outside the forward shock front at a low 
speed of ~ 10 km s~^. As in LZ , we adopt spherically sym- 
metric nonlinear hydrodynamic partial differential equations 
(PDEs) in spherical polar coordinates (r, 9, ip) to model the 
piecewise isothermal gas dynamics of these four wind regions 
of a quasi-spherically symmetric PN. 



2.1 Inner Shocked Tenuous Wind Zone 
Featuring a Reverse Shock 

As already noted, a reverse shock emerges and propagates 
in the inner fast wind zone due to the coUisional interaction 
between the later fast and earlier slow winds. The mass and 
radial momentum conservations across an isothermal shock 
front in the shock comoving framework are simply 

pd{Ud — Us) = pu{Uu — Us) , 
<hiPd + PdUd{Ud — Us) = a^pu -I- puUu{Uu — Us) , (1) 

where a is the isothermal sound speed, p is the mass density 
and u is the radial flow velocity, and the two subscripts d and 
u denote the downstream and upstream sides of an isother- 
mal shock an d the subscript s indicate association with 
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the shock (e. g. ICourant fc Friedrichs l 1976; Spitzer 
Shen fc Lou I l2004l : iBian fc Lou I . |2005; Yu et all 
LZh . The temperature ratio across s uch a shock front i s 
another parameter to be specified fsee lShen fc Loulliooi ). 
The reverse shock travels inwards relative to the inner wind, 
so the inner fast wind region is the upstream side of the 
shock and the reverse shock heated wind occupies the down- 
stream side until the contact discontinuity at the expand- 
ing radius R^. In the laboratory framework of reference, 
a reverse shock may move inwards or outwards, or appear 
stati onary in space dep ending on various specific situations 
(e.g.,[LE3,[l99i,[l99i). 

We presume a constant mass loss rate Mf^ for an 
isothermal inner fast wind and define a^^ci(u) , T^^d{u) , Pw,d{u) 
and as the isothermal sound speed, constant tem- 

perature, mass density and constant wind velocity on the 
downstream (upstream) side of the reverse shock. The sound 
speed ratio r„ = a,„_d/a„,„ = {T^,d/Tw,u)^^'^ across the 
shock front characterizes the downstream plasma heating of 



the reverse shock. Isothermal shock conditions ^ then give 

Vw,u — Vrs = {l/2)[Vnj.d ~ Vrs + O^w .d / {'"■m .d — Vrs)] 



+ 



I \2 2 

\V-m,d-Vrs) - au,,d 

{Vw,d — Vrs) 



1/2 



+ iaw,d — 



Vw,d — Vrs = (l/2)[Uu,,u — Vrs + ,ii / (^lu ,u — Vrs)] 



{Vn,,u - Vrs) ' 



2 -,2 



1/2 



Pw,d{RrB,t) = Pw,u{Rrs,'t){Vw,u ~ Vrs)/{Vw,d — Vrs) , (2) 

where Vrs is the rev erse shock velocity in the laboratory 
reference framework ||LZ I) . Calculations show that positive, 
negative and zero Vrs are all physically allowed. For Vrs < 0, 
the inner fast wind region shrinks radially inwards until a 
certain epoch. In other words, the reverse shock heats up the 
downstream wind plasma and leads to a hot bubble inside 
the contact discontinuity surface at an expanding radius Re- 
The mass density profile within the central fast wind 
zone is simply 



pm,u{r, t) = M/„/(47rw„,ur^) 



(3) 



and that in the reverse-shocked downstream wind zone is 



Pw,d{r,t) 



1 — Vrs/Vn,,u Mfw 

{vw,d - Vrs) 47rr2 



Mf^r- 



47r(«»,d 



(4) 



which remains steady, where the constant (3 parameter is 
explicitly defined. This last approximation highlights that 
the downstream density of shocked inner wind does not sen- 
sitively depend on the inner fast wind speed Vw,u because 
we typically have a small speed ratio of Vrs/vw,u ~ lO"'^. 
Due to the reverse shock, the downstream mass loss rate 
just inside the contact discontinuity surface is given by 

M' = 47vR^,v^,,p^ARc, t) = ii^^^!^M^„ , (5) 

(1 — Vrs/Vru,d) 

which is equal to the central mass loss rate M only when 
the reverse shock remains stationary in the laboratory refer- 
ence framework. The difference between Mf^ and M' is the 
mass loss rate at the reverse shock front due to the reverse 
shock movement. 

X-ray observations in the past several years appear to 
indicate that the temperatures of inner shocked fast wind, 
where X-ray photons are released, are a factor of ~ 2 — 10 
lower than that predicted by simple energy conservation 
arguments (e.g. Soker & Kastner 2003). In our formula- 
tion, we use a chosen sound speed ratio r^, = a^^d/o.w,u = 
(T^.d/Tui.u)^^^ to bridge the two constant temperatures of 
upstream and downstream sides across a moving reverse 
shock. 



2.2 The Dense Shell Region and the 
Outer AGB Slow Wind Envelope 

The shocked dense shell region and the outer AGB slow 
wind envelope separated by an outgoing forward shock are 
described by the ISSV model of ,LZ . . The piecewise isother- 
mal hydrodynamic equations with spherical symmetry are 
the mass conservation 

^^=4.r'p, (6) 



dM dM „ 



dr 
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or equivalently, 

dp I d . 2 N 
Ot or 

and the radial momentum equation 



du du 
dt dr 



a^dp 
p dr 



= 



GM 



(7) 



(8) 



where u is the radial bulk flow speed; M(r, t) is the enclosed 
mass within radius r at time t; p{r, t) is the mass density; 
a = [p/pY^"^ = [k-oT /ffiY^'^ is the isothermal sound speed; 
T is the constant gas temperature; m is the mean particle 
mass; p is the gas pressure, fee is Boltzmann's constant and 
G = 6.67 X 10~* cm"^/(g s^) is the gravitational constant. 
By the known isothermal self-similar transformation 



X = r /{at) , 
M{r, t) = a^tm{x)/G , 
equation ^ leads to 

m(x) = x^a(x — v) 



pir, t) = a{x)/(4nGt^) , 
u{r, t) — av{x) , (9) 



(10) 



and then PDEs ((?]) and ([8]) reduce to two coupled nonlinear 
ordinary different equations (ODEs) 



[{x- 



1 da 
a dx 



a{x - 



{x — v) 



(x — v) 



(11) 



(12) 



where a{x), m(x), v(x) are the dimensionless reduced 
variables corresponding to mass density p(r, t), enclosed 
mass M{r, t) and radial flow speed u(r, t), respe c tively 
llShu I. Il977l: iHunter I. Il977l: IWhitworth fc SummeiTI. |l985: 



Tsai fc Hsu] 



Shen fc Lou I . 



1995 



2004 



Shu et al. 



Bian fc Lou I 



20021: Lou fc Shen I , 12004 : 
2003 : |LZ| . |2009| ). Given 



proper analytic asymptotic solution conditions at large and 
small X as well as in the neighborhood of the sonic critical 
point, these two coupled nonlinear ODEs pi|) and H12|) can 
be solved in a straightforward manner. Each dimensionless 
ISSV solution of these two ODEs gives a behaviour of the 
dense shell and outer AGB wind envelope outside Rc- 

The forward shock travels outwards in a self-similar 
manner and divides the shell into the shocked dense shell 
region and the outer slow AGB wind envelope. We define 
ad(u), Ta(u), ad(u) and Va(u) as the sound speed, tempera- 
ture, reduced mass density and reduced gas radial flow ve- 
locity on the downstream (upstream) side of the forward 
shock and the sound speed ratio t = ad / au = (Td/Tu)^^^ 
JShen fc Lou I. l2004l : iBian fc Loul l2005l : IYu et al. I . |2006| . 
ILZ h . In this case, the dense shell region is the downstream 
side and the outer slow AGB wind envelope is the upstream 
side of the outgoing forward shock front. We define dimen- 
sionless variables Xsd = Vs/ad and Xsu = Vs/a-u where 
Vs = Tsjt is the radially outgoing velocity of the forward 
shock in the laboratory reference framework. Then isother- 
mal shock conditions ([l]) bear the self-similar form of 

a<i/Q« = (wu - Xsv)l{T(vd - Xsd)\ , 

(Vd - Xsd)/{Vu - Xsu) - T = {rVd - Vu){Vd - Xsd) (13) 

(Shen fc Lou 2004). These isothermal shock jump conditions 



connect the AGB outer wind envelope with the dense shell 
region across the forward shock front. 

The analytic asymptotic solution of coupled nonlinear 
ODEs (lllf) and (|12|l at large x is known as 

^._y I 2-A ^ V ^ {A/6-l){A-2) + 2V^/3 ^ ^ 
x x'^ X? ' 

A A{2-A) U~A)VA 

where V and A ar e the velocity and mass parameters 
(|Lou fc Shen I |2004| : ILZ h . In our self-similar piecewise 
isothermal model framework, these two parameters charac- 
terize the asymptotic massive AGB wind velocity as 

t^AGB = Va„ (15) 

at very large x and a constant AGB wind mass loss rate as 



Magb = VAallG . 



(16) 



These asymptotic relations are used to constrain the range 
of model parameters for PNe. 

2.3 Radial Expansion of the Spherical 
Contact Discontinuity Interface 

The contact discontinuity surface at Rc expands at a con- 
stant speed and separates the inner and outer zones. ILZ I 
estimated that the gravity of inner tenuous wind zone may 
be negligible during the dynamic self-similar evolution of 
the dense shell region and outer AGB envelope in outer 
zone. By letting m{xo) = in equation (llOp . we then ob- 
tain the contact discontinuity surface radius Rc — adXot and 
its expansion velocity Vc = UdXo- By doing this, we ignore 
the self-gravity from the inner shocked wind region within 
the contact discontinuity radius Rc, but still retain the self- 
gravity of the dense shell and the AGB slow wind envelope 
outside the contact discontinuity radius Rc- A contact dis- 
continuity requires that the inner downstream wind velocity 
Vuj,d and the outer dense shell gas velocity at Rc be equal to 
the contact discontinuity expansion velocity, namely 

Vuj,d = adXo = adv{xo) , (17) 

which is satisfied automatically when m{xo) = occurs. 
Here, xo is where an ISSV solution for the dense shell re- 
gion and the outer envelope starts. Given a proper set of 
{a;o, ao = a{xo), Xds}, one can readily integrate coupled 
nonlinear ODEs (|ll|l and (|12|l from xq to larger x > xo (i.e. 
from left to right), apply self-similar shock jump conditions 
([13} at Xds to cross an outgoing forward shock front once 
ad and au are known and then continue further from Xus to 
sufficiently large x to determine velocity and mass parame- 
ters V and A with desired accuracies. The flow zone between 
the contact discontinuity and the forward shock is the dense 
downstream side and the outer AGB slow wind is the up- 
stream side. The forward shock jump condition connects Xds 
and Xus- These solutions have a negligi ble m ass within xo 
and are referred to as ISSV solutions bv lLZ L In total, eight 
parameters are involved for constructing a model ISSV so- 
lution {xo, ao,ad, Xds, au,Xus, V, A}, with five of them being 
actually independent. 

At radius Rc of contact discontinuity, the inner down- 
stream wind plasma pressure is 



Pw,d{Rc, t) 



kBTui,d P ^ Tu),d 

rfiuiM Rl 



(18) 
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where fh^^d is the mean particle mass of the inner down- 
stream wind zone. With self-similar transformation ((9|, we 
have the shocked dense shell gas pressure at the contact dis- 
continuity radius Rc as 

^''(^-') = ^4^°^F' ^^^^ 

where rfid is the mean particle mass of the shocked dense 
shell zone. To physically maintain an expanding contact dis- 
continuity at constant speed, the parameters M/^, Vrs, 

ao and Td are coupled by the pressure balance 
condition Pw,diRc,t) = Pd{Rc,t) across the contact discon- 
tinuity, namely 

Td ap _ T^,d 47r/3 _ Tn,,d (1 — Vrs/vw,u) Mf^ 
fhd G fh^^d v^^^ fh^^d (1 - Vrs/v^^d) 

In general, the two mean particle masses fhd and m„,tj are 
allowed to be different depending on specific situations. 



2.4 Comparison with Previous Models of PNe 

IChevalier fc Imamura I ()l983h and ILZ I both invoke a self- 
similar hydrodynamic wind-interaction phase to explore the 
global dynamic shock evolution of PNe. In both models, the 
inner and outer wind zones are characterized by different yet 
constant mass loss rates. 

More specifically, in the model of Chevalier & Ima- 
mura (1983), the innermost fast wind zone and the out- 
ermost AGB slow wind envelope are both isothermal and 
expand at constant speeds with both mass densities scaling 
as oc . In between these two wind regions, both shocked 
downstream wind zone and dense shell are solved separately 
by adopting a type of adiabatic self-similar transformation 
for 7 = 5/3 hydrodynamic PDEs without gas self-gravity 
(Parker 1961). These two self-similar dynamic zones are 
joined by the pressure balance across an outgoing contact 
discontinuity. Two different Mach numbers are specified re- 
spectively for the reverse shock separating inner fast 
wind and self-similar shocked downstream wind zone (hot 
bubble), and for the forward shock separating self-similar 
dense shell and out er slow AGB wind envelope. 

In comparison, ILZ I model also presumes a density 
profile for a steady piecewise isothermal fast wind zone in 
the central region around the remnant white dwarf. Dif- 
ferently, the dense shell and the outer AGB slow wind 
envelope are described by a self-similar dynamic solution 
with self-gravity and an is other mal shock jump (i.e. type 
Z self-similar solutions of ILZ I ). The self-similar process 
here differs from that of Chevalier & Imamura (1983) in 
(a) the inclusion of self-gravity; (b) the piecewise isother- 
mal wind (both the dense shell and the outer AGB slow 
wind envelope are isothermal but with different tempera- 
tures); (c) a different self-similar transformation due to the 
presence of self-gravity and (d) necessary shock parame- 
ters (forward shock is described by its velocity and sound 
speed ratio or temperature jump param eter r). In the in- 
ner shocked downstream wind zone. IlzI model has constant 
temperature and wind speed. For a constant mass loss rate 
M/u, = Pu,,dWm,d47rr^(l-Wrs/u»,d)/(l-Wri,/fi£.,u), both mass 
density and thermal gas pressure scale as r~'^ . 

In the inner shocked downstream wind zone and the 



dense shell, ILZ I piecewise isothermal model behaves dif- 
ferently from adiabatic self-similar solutions of Chevalier 
& Imamura (1983). For Chevalier & Imamura solution of 
colliding winds, mass density and pressure in both down- 
stream zones (i.e., the inner shocked downstream wind and 
the dense shell) are lowest near the shock fronts and high- 
est near the contact discontinuity where strong shocks were 
initiated. Under weak shock condition, density behaves sim- 
ilarly but pressure is no longer mo notonic, which may lead 
to R ayleigh- Taylor instability (e.g. IChevalier fc Imamura I , 
Il983l ). Several isothermal cases with 7=1 for the dense 
shell are also calculated by Chevalier & Imamura (1983) 
to find that both density and pressure peak at the forward 
shock front, instead of at th e contact discontinuity, which is 
qualitatively similar to lLZ I solutions for the dense shell. An- 
other difference is the behaviour of AGB envelopes. In ILZ I 
model, both velocity and mass loss rate of AGB envelopes 
approach constant at large radii. 

Based on Chevalier & Imamura (1983) adiabatic solu- 
tion for inner shocked fast wind, Akashi et al. (2006) find 
the morphology of the X-ray emission is in the form of a 
narrow ring just inside the optical bright part of a PN (i.e. 
dense shell). However, simulation shows that the velocity 
of reverse shock and contact discontinuity cannot stably re- 
main constant (Akashi et al. 2007). This might be related 
to the Rayleigh- Taylor instability for some adiabatic solu- 
tions of Chevalier fc Imamura (1983). Akashi et al. (2007) 
also find that in order to get more realistic simulation re- 
sults, a rapidly decreasing mass loss rate of inner fast wind 
is required. In this case, the density should be no longer 
increasing from the reverse shock to the contact discontinu- 
ity. For Chevalier fc Imamura isothermal solutions and for 
ILZ I model, mass densities increase monotonically from the 
contact discontinuity to the shock front. A physical conse- 
quence is that near the reverse shock gas is brighter in X-ray 
emissions than near the contact discontinuity. 

For a PN, both forward and re verse shocks are radiative 
(e.g. IChevalier fc Imamura I . 1 19831 ). The central hot white 
dwarf can be a strong source of photoionizing radiation 
and may compl etely photoionize the outer dense shell (e.g. 
Chevalier 1997; ILZ h . Therefore, in reality, the hot bubble 
and dense shell should involve gas dynamics between adia- 
batic and isothermal processes (say, a polytropic description 
as an approximation). 



3 DIFFUSE X-RAY EMISSIONS FROM PNE 

The reverse shock initiated by the collision of an inner fast 
tenuous wind with ambient dense gas materials previously 
blown out during the AGB slow wind phase heats up the 
inner downstream wind region, leading to spectral line emis- 
sions from highly ionized metals such as C, O, Ne in X-ray 
bands (e.g. Volk fc Kwok 1985; Akashi et al. 2006; Leahy 
et al. 2000; Stute fc Sahai 2006; Yu et al. 2009). In con- 
trast, temperatures of the shocked dense shell and of the 
outer AGB wind envelope are onl y in the order of ~ 10* 
K (e.g. ICuerrero et al. 1. 120041 : ILZ h . and that of the central 
hot white dwarf is < 2 x 10^ K (e.g. Werner et al. 2008). 
Consequently, X-ray emissions from PNe would be mainly 
produced in the inner shocked downstream wind region with 
a higher temperature. In this scenario. X-ray emissions from 



6 Y.-Q. Lou & X. Zhai 



PNe projected in the plane of sky should largely lie within 
the expanding co ntact discontinuity of radius Rc- 

According to ILZ I . the typical electron number density 
in the dense shell zone is estimated as rie ~ 10^(1000yr/t)2 
cm~'^, where t is the typical age of a PN, and the thickness 
of a dense shell zone is estimated as ~ 10^*(f/1000yr) cm. 
Then the mean free path of a photon in the dense shell zone 
would be ; = l/(neO-) ~ 10^''(t/1000yr)2 cm, where a = 
6.65 X 10~^* cm^ is the electron cross section for Thomson 
scattering. Therefore the dense gas shell should be optically 
thin for emissions of X-ray photons during the evolution of 
a PN. We denote the X-ray emissivity of the shock heated 
downstream inner wind plasma by A(r) such that Anp is 
the energy emitted per unit time per unit volume, where 
rip is the proton number density (e.g. Sarazin 1986). At the 
temperature range of T < 3 x 10^ K when spectral line 
emissions from highly ionized metals are important, a simple 
approximation is thus A(r) « 6.2 x IQ^T"" '' ergs cm"^ 
sec~^ (e.g. Sarazin 1986). The X-ray luminosity Lx of a PN 
can then be calculated by integrating the emissivity over the 
entire shocked downstream inner wind region, viz. 



Lx = 



(21) 



where is the proportion of gas mass density to proton 
number density. For a H II region, we simply have m,i = nip, 
while rUfi = 1.18mp for gas with a helium abundance of 
He/H= 0.1/0.9 by number. Defining the X-ray parameter 
X = A/3^/m^, we readily obtain the following expression 



Lx =47rX(l/i?,, - l/Rc) 



(22) 



A normalized X-ray surface brightness profile B can be cal- 
culated by integrating over the X-ray luminosity along the 
line of sight. Let R be the projected radius from the centre of 
a grossly spherical PN, we have the X-ray surface brightness 
radial profile B as 



A 



dx , 0<R<Rr 



dx , Rra < R < Rc 

JO rn^ 

. R> Rc 

(23) 

By introducing a dimensionless function L{z) defined as 

L{z) = [^(1 - z^f^'^ + arccosz] for Q<z<l , (24) 
we immediately arrive at the following expression 

= < . I R 



B{R,t) 
A/i?3 



L 



, 



Rc 



Rrs ^ R K Rc 

R> Rc 



(25) 

Figures [T] and [2] illustrate several model PN X-ray bright- 
ness B profiles marked with different values of radius ratio 
Rrs/Rc within the range for this normalized reverse shock 
radius < Rrs/Rc < 1. We emphasize that morphologies 
of such diffuse X-ray emissions from PNe are determined 
by this radius ratio Rrs/Rc- The case of Rrs/Rc close to 
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Figure 1. Shown here are variations of the projected dimension- 
less X-ray brightness radial profile i3/[A/3^/(m^i?^)] within the 
contact discontinuity radius Rc for PNe with different radius ra- 
tios Rrs/Rc- The projected radius R is also normaUzed by the 
contact discontinuity radius Rc- The range of brightness B vari- 
ation spans three order of magnitudes as Rrs/Rc varies from 0.1 
to 0.9. Different epochs of evolution might be invoked to explain 
why some PNe are observed to be strong X-ray sources while 
other sources of PNe are not. 



unity corresponds to a ring-like brightness profile (e.g. the 
case of PN NGC 40), while a small Rrs/Rc gives a more 
spherical appearance for the X-ray surface brightness radial 
profile with a more luminous central bulge or sphere (e.g. 
the case of PN BD -f 30° 3639). Both X-ray luminosity and 
surface brightness radial profile depend only on the X-ray 
parameter X and radius ratio for reverse shock and con- 
tact discontinuity. The X-ray parameter X is a function of 
A (or inner shocked downstream wind plasma temperature 
Tw,d), P and while parameter jS is further determined 
by the mass loss rate Mf w from the central white dwarf, 
reverse shock speed Urs, inner fast wind velocity Vw,u and 
shocked downstream wind velocity Vw,d according to expres- 
sion (Q. As the inner radius of the dense shell, the contact 
discontinuity radius Rc can be measured using observations 
in optical and infrared bands (e.g. Kastner et al 2000; Guer- 
rero et al 2004; Montez et al. 2005). Parameters X and Rrs 
can be inferred by fitting an X-ray surface brightness radial 
profile and the estimate of total X-ray luminosity of a PN. 
Spectral analysis can measure the innermost fast wind ve- 
locity, dense shell gas velocitjQ and temperatures (four in 
all) in the inner and outer zones (e.g. Guerrero et al 2004; 
Montez et al. 2005). Pressure balance condition (|20[) across 
the expanding contact discontinuity radius Rc further con- 
nects physical parameters in the inner shocked downstream 
wind region and the dense shell. We emphasize that com- 
bined observations in infrared, optical and X-ray bands may 
be utilized to infer all relevant parameters that characterize 
a quasi-spherical PN. 

The morphology of diffuse X-ray emission from a PN 

^ Relation (I17I I indicates that the dense shell gas velocity also 
gives the inner shocked downstream wind velocity as required by 
the pressure balance across an expanding contact discontinuity 
radius Rc- 
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Figure 2. Gray-scale plots of the logarithm of dimensionless X- 
ray surface brightness (see Fig. [TJ with different radius ratios 
0.1 < Rrs/Rc < 0.9 marked at the upper-right corner of each 
panel (see Fig. [TJ. The abscissa and the ordinate represent nor- 
malized projected radius R/Rc by the contact discontinuity ra- 
dius Re- As the radius ratio Rrs/Rc may vary with time, the 
corresponding variations in brightness profile is thus expected. 



only depends on the radius ratio Rrs/Rc (see Figs. 1 and 2). 
Thus the time evolution of this radius ratio determines how 
the X-ray morphology evolves with time t. In other words, 
different epochs of PN evolution may be invoked to explain 
why some PNe are observed to be strong X-ray sources while 
others are not. Here, we present straightforward calculations 
to show the relation between the physical evolution state of 
PNe with their diffuse X-ray emission morphologies. Sup- 
pose the contact discontinuity surface has a radius Rco at 
the time when the central fast wind collides with the slow 
dense AGB wind envelope. Let t* be the time lapse from 
when the time forward and reverse shock pair initiated after 
a wind-wind collision to the current epoch. We may use Rco 
and Vc to estimate a kinematic timescale t ~ Rco /vc + t* 
and regard this timescale t as the dynamic age of a PN (see 
self-similar transformation equation[9]above). It follows that 
it — f) roughly represents the temporal duration that the 
central star has depleted its hydrogen envelope but not yet 
started to launch a fast hot tenuous wind. Then the radius 
ratio Rrs/Rc evolves with time as 

Rrs _ [RcO + Vrst*) _ Vet ~ {Vc ~ Vrs)t* 



Vet 



(26) 



[RcO + Vet*) 

At the moment that wind-wind collision occurs and both 
forward and reverse shocks are initiated, the radius ratio 
Rrs/Rc = 1 and no diffuse X-ray is emitted (see egn. I22p . 
If PNe are identified around this stage, then there would 



be no or very weak X-ray emissions. The time evolution of 
Rrs/Rc is mainly determined by Vrs according to expression 
(I26|) . If Vrs is nonnegative, then this radius ratio has a lim- 
iting value of Vrs/vc in the long term evolution for a large 
t* (see Fig. 3) during which the initial contact discontinu- 
ity radius Rco becomes less important. Otherwise, when the 
reverse shock travels radially inwards, the region inside the 
contact discontinuity surface tends to be filled with more 
shock heated downstream wind plasma. Meanwhile, grav- 
ity, magnetic field and radiation of the central star at the 
shrinking radius of the reverse shock will eventually become 
sufficiently strong, and neither asymptotic central wind pro- 
file nor shock condition ([1} is reliable any more. We will 
show presently that the energy conservation actually forbids 

Rrs 0. 

We now consider this problem from the perspective of 
energy conservation. The inner fast wind kinetic energy and 
thermal energy input per unit time 

dEin 1 2 I 3 • 2 /o'7^ 

= -^Mf^v^^^ + -Mf^a^^^ (27) 

should be larger than the kinetic and internal energy in- 
crease of the system when X-ray is emitted. The internal 
energy change rates of fast wind zone and shocked fast wind 
zone are 



dt 



,u 3 ,^2 /'r>^2 3 Vrs , , 
— — -4:TrRrsVrsPw,u(Rrs)aiu,u ~ TT -M f 



2 Vw,u 



dt 2'^'^RcVcPw,d{Rc)a^^a - -'inRrsVrspw,d{Rrs)a^^eL 



= -(1 — Vrs/Vw,u)MfwChn,d 



(28) 



respectively. The kinetic energy change rate caused by the 
reverse shock in the inner wind is 



dK 
~dt 



AtvR; 



iVw,d—Vrs) 



Pw^dV^^d 



'iVw,u Vrs) 



/I / N / 2 2 X , n 

——(1 - Vrs/Vyj,u)(v^^a ~ < 



(29) 



The power of the shocked downstream wind zone (hot bub- 
ble) pushing the dense shell outwards is 

dW 2 2 ,> 2 (l-Vrs/v^^^) 

—— — 'nV-tteV^ dPw.dO-m.d — Mf^a^ j j- ; r- . 

at [I — Vrs/V^^d) 

(30) 

Given typically estimated parameter values of a PN 
Vw,u ~ 10^ km s~^, Vya,d = Vc ~ 10^'"^ km , Vrs ~ lO""'^ 
km s~^, ai a.u ~ 10 km s~^ and a^^d ~ 10^^^ km s^^ as 
inferred bv lLZ I . we find that 

dEin ^ ^q2^3 dUw^d dUio,d ^ dUw,u 



dt dt 
dEin dK 
~~dt df ' 



dt 

dEin _ 10^-^3'^ 
dt ^ dt 



dt 



(31) 



The energy conservation then requires that 

Lx < ^ - ^ {U^,u + U^,d + K + W)^ Mfn.vl,u- (32) 

Energy conservation requirement 1)32^ must be neces- 
sarily met. Using expression ((22]), we find that equation 
(|32[) limits the possible range of the hot bubble tempera- 
ture T^i.d and the reverse shock radius Rrs- As the X-ray 
luminosity Lx approaches the upper limit M/u,ii^_„, either 
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Figure 4. Multipole expansion to X-ray counts distribution of 
NGC 40. Equivalently, this multipole spectrum is the Fourier 
transformation of the polar angle distribution of X-ray brightness 
map into </>— space, where <j) S [— tt, -|-7r]. The x— axis is the index 
of multipole and the y— axis is the normalized absolute value of 
multipole spectrum so that the central strength is 1. 



Figure 3. Temporal evolutions of the radius ratio Rrs/Rc (Panel 
A), the total X-ray luminosity (Panel B) and the average X-ray 
surface brightness (Panel C) of three illustrating tracks for PNe 
with different reverse shock velocities Vrs- The solid, dotted and 
dash curves in three panels show how PNe evolve when Vrs = —3, 
and 10 km s~^, respectively. All three models have the same 
innermost fast wind mass loss rate Mfw = 1O~*M0 yr~^, contact 
discontinuity radius expansion speed Rc = Vc = 30 km s~^, inner 
shocked downstream wind plasma X-ray emissivity A = 10"^^ 
ergs cm^ and R^o = 1-3 X 10^® cm (the distance that a wind 
blows by at 20 km s~^ for ~ 200 yrs). An inner fast wind velocity 
of ~ 500 km is adopted here. 



the cooling becomes dominantly important so that the tem- 
perature r,„_d drops rapidly, or the reverse shock radius Rrs 
can no longer shrink. However, in reality, it is unlikely that 
X-ray emissions from the hot bubble reach the upper limit 
because cooling efTects would become extremely important 
before Lx ~ M/mii^ In the above discussion, we have ig- 
nored the gravitational potential energy change rate of the 
shocked downstream zone occupied by the hot bubble. 

Fig. [3] illustrates how the radius ratio Rrs/Rc, the to- 
tal X-ray luminosity and the average surface brightness of 
PNe evolve with negative, zero and positive reverse shock 
velocities Vrs- While all three conditions give decreasing ra- 
dius ratio Rrs/Rc in general, the temporal evolution of total 
X-ray luminosity and average surface brightness differ from 
each other significantly. In the case of Vrs ~ —3 km s^^, 
the reverse shock radius Rrs stops at ~ 10 an as required 
by energy conservation (|32p . We will show presently that for 
both PNe NGC 40 and NGC 7662, X-ray emissions only take 
up < 1% of the inner fast wind energy input. Calculations 
show Rrs ~ 1000 an when Lx reaches ~ 1% of M/^n^ 
viz. Lx = 10'^^ erg s~^. This appears consistent with the 
earlier estimate of Frankowski & Soker (2009). 



4 OBSERVATIONS OF SAMPLE PNE 
4.1 Planetary Nebula NGC 40 

NGC 40 (PNG120.0+09.8; lAcker et al.1ll992l ) is a very low 
excitation PN powered by a WC8 star evolving towards a 
white dwarf. Within a ring of nebulosity revealed by optical 
and near-infrared images of NGC 40, Montez et al. (2005) 
detected faint, diffuse X-ray emission distributed as a partial 
annulus using the data from the Chandra X-Ray Observatory 
(ObsID 4480). This annulus-like X-ray profile corresponds to 
a fairly large ratio of Rrs/Rc and thus a low X-ray luminos- 
ity (see equations 1221 and 1251 and Figs. [T] and [2]). The actual 
inferences are X-ray temperature T = (8.0 ± 2.0) x 10''' K 
and total luminosity Lx ~ 1.5 x 10^° (D/l.Okpc)^ ergs s~^ 
of NGC 40 (Montez et al. 2005). Kastner et al. (2008) repro- 
cessed the X-ray data of NGC 40 and updated the luminos- 
ity value to Lx ~ 4.0 x lO^'^erg s"'^ which is adopted in our 
model calculations. The X-ray temperature and luminosity 
are one of the lowest measured so far for PNe with diffuse 
X-ray emissions. Here, the distance to NGC 40 -D = 1.0 kpc 
is estimated by Leuenhagen et al. (1996). 

We decompose the X-ray counts distribution in terms 
of cylindrical harmonics in polar coordinates in the plane 
of sky and demonstrate the resulting multipole spectrum in 
Figure|4]above. This spectrum reveals a strong polar uniform 
component (multipole index — 0) and a distinct dipole com- 
ponent (multipole index — 2) which is about one third times 
as strong as the uniform component. While being somewhat 
crude, this suggests an exploration for the X-ray brightness 
radial profile of NGC 40. 

The observed X-ray brightness radial profile of NGC 40 
is displayed in Fig. [5] (i.e. the open circles with error bars). 
We assume a uniform X-ray emission background approxi- 
mately and fit our PN model results (i.e. the dotted curve 
in Fig. [5]) with the observed data profile. The best-fit model 
gives a radius ratio of Rrs/Rc ~ 0.8 and an X-ray parame- 
ter X = 2.83 X lO''* {D/1.0 kpcf ergs s"^ cm. By equation 
(|22|l . these results further lead to an X-ray luminosity from 
NGC 40 as Lx ~ 3.7 x 10^^(I>/1.0 kpc)^ ergs s'S consis- 
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tent very well with the observational inference (e.g. Kastner 
et al. 2008). The model fit also gives a reverse shock radius 
Rrs « U"D = 2.0 X 10"(I>/f.O kpc) cm and a contact 
discontinuity radius Rc ~ 17" D = 2.5 x fO^^(D/f.O kpc) 
cm, corresponding to the peak and boundary of diffuse X- 
ray emissions from shocked hot bubble. In Figure [H we cir- 
cle these model-fitted Rra and Rc on the optical and X-ray 
composite image of PN NGC 40. By visual inspection, Rr3 
does pass through the X-ray brightest region and - Re ro ughly 
marks the boundary of the diffuse X-ray emission. ILZ I piece- 
wise isothermal model claims that the optically bright PN 
rim is located just outside the contact discontinuity radius 
Rc and the density becomes highest near the forward shock 
radius Rs . Figure [6] indicates that the forward shock is just 
outside the contact discontinuity surface so that the optical 
bright dense shell is very thin, i.e. -Rs > -Rc ~ 17". Accord- 
ing to the spectral analysis of Sabbadin et al. (2000), the 
most optically luminous shell region of NGC 40 still have 
redshift and blueshift componeonts at radius > 17", imply- 
ing Rs > 17". Figure[5]shows that our model fits observation 
very well. The radial separation between the most optically 
bright radius at ~ 17" and the most X-ray bright ra dius 
at ~ 14" is clearly distinguishable, strongly supporting ILZ I 
piecewise isothermal model. In contrast, the adiabatic model 
of Chevalier & Imamura (1983) would claim no significant 
separation between these two radii (e.g. Akashi et al. 2006). 

The best-fit Raymond-Smith thermal plasma emission 
model (Raymond & Smith 1977) with the X-ray spectrum 
of NGC 40 indicates a plasma temperature of Tx = (8.0 ± 
2.0) X lO'^ K (Montez et al. 2005); at such temperature, the 
spectral line cooling from highly ionized metals is impor- 
tant while the thermal bremsstrahlung emission from free 
hot electrons is relatively weak, and the total X-ray emis- 
sivity is A « 2.0 x 10~^^ ergs cm^ s"-^ (e.g., Sarazin 1986). 
With rUfj, = mp, we have P = 3.16 x lO"*M0yr" V(100km 
s~^) by expression Q. We adopt Vw,d = Vc — 30 km s^^ 
as estimated by optical spectrum line analysis near the con- 
tact discontinuity surface (e.g., Sabbadin et al. 1999). Ignor- 
ing the duration that the central star stopped to blow slow 
AGB wind and not yet began to launch tenuous fast wind, 
we estimate a wind-wind dynamic interaction timescale of 
Rc/vc ~ 3000 yrs (a comparable timescale of ~ 5000 yrs 
was estimated by Akashi et al. 2006). For the radius ra- 
tio Rrs/Rc = 0.8, we have Vc — Vrs = 6.0 km s~^. Thus 
the reverse shock moves outwards in the laboratory frame- 
work of reference. Meanwhile, equation (2) gives an inner 
tenuous fast wind speed v^.u = 1.1 x 10'^ km s~^. Equa- 
tion Q gives a central star mass loss rate as Mf^ ~ 
2.4 X 10~* M0yr~^, which is very close to an earlier esti- 
mate of ~ 3 X 10~*Mq yr~^ by UV spectral analysis, and 
falls within the estimated range of ~ 10~^° - IQ-'^Mq yr-i 
(e.g., Cerruti-Sola & Perinotto 1985). By imposing condi- 
tion (|20p . the diffuse X-ray emission from the central zone 
of PN NGC 40 requires a Td = 1.5 x 10''(10"Vao) K for the 
outer shocked dense shell. We note that the X-ray emission 
only takes up ~ 0.2% of M/™-u^,„ for PN NGC 40. 

However, our best-fit model in Fig. [S] overestimates the 
X-ray brightness for smaller radii around the center. In lLZ I 
model, the temperature of the dense shell region is not high 
enough to radiate in X-ray bands. In reality, thermal con- 
duction inevitably occurs across the contact discontinuity so 
that the dense shell region will be heated up by inner shocked 



downstream wind and may lead to partial X-ray emissions. 
The inner shocked downstream wind zone near the contact 
discontinuity surface would also become denser and cooler 
due to diffused materials and thermal conduction from the 
dense shell, especially if magnetic field is absent (e.g., Steffen 
et al. 2008). This effect would lead to a relatively higher X- 
ray flux than our model result around the outer boundary of 
the hot bubble. Nevertheless, this thermal conduction is ex- 
pected to be sensitively suppressed even for the presence of 
a weak magnetic field transverse to the radial direction (e.g.. 
Chevalier 1997; Yu et al. 2006; Wang & Lou 2008; Lou & 
Hu 2009: [lzI ). so the total X-ray luminosity of PNe is more 
or less confined to emissions within radius Rc- Within the 
radius of ~ 2 arcsec or ~ 3.0 x 10^^ cm or ~ 2000 au around 
the central region. X-ray brightness of our model is ~ 75% 
greater than the upper limit of confidence interval of obser- 
vation (see Fig. [5]). Since the reverse shock is at ~ 14" within 
which there is no diffused X-ray emission, the projected X- 
ray photons within radius ~ 2" are actually emitted in the 
hot bubble at larger radii closed to Rc but projected around 
the center along the line of sight. If there is a dense dust 
cloud surrounding the central star, the X-ray photons emit- 
ted behind the cloud might be partially absorbed so that 
the observed central X-ray brightness becomes lower than 
expected. Both theoretical calculations (e.g., Frankowski & 
Soker 2009) and observations show that dense gas/cool dust 
clouds may exist around the central white dwarf of PNe with 
diameter > 150 au (e.g. Helix nebula, dust density ~ 2.5 
g cm~'^ by Su et al. 2007). The cross section of such a dust 
cloud is about 2 x 10* au^, which can only obscure ~ 0.5% 
X-ray emission within ~ 2000 au radius maximally. If this 
is the case for PN NGC 40, Figure [5] requires a column elec- 
tron number density of ~ 2.7 x 10^* cm~^ within ~ 2" to 
absorb ~ 80% X-ray emissions behind central star by Thom- 
son scattering. Such a high column density corresponds to 
a total cloud mass ~ 5Mq. We thus conclude that X-ray 
absorption by either dust or gas cloud cannot explain the 
unusually low X-ray brightness in the central NGC 40. Al- 
ternatively, the non-spherical structure of NGC 40 might be 
the most important factor. 

The dipolar structure of the X-ray brightness is no- 
ticeable in PN NGC 40. Also, the eastern X-ray region 
is brighter and seems to be more confined to within the 
brighter optical rim. As a first approximation, our spher- 
ically symmetric treatment modelling here ignore these 
asymmetries. 



4.2 Planetary Nebula NGC 7662 

NGC 7662 appears to be a quasi-spherical PN with a mod- 
erately large diameter of ~ 30" or ~ 4 x 10^^ cm and with 
optical surface brightness but without X-ray data available 
up to the present. Guerrero et al. (2004) probed the struc- 
ture and kinematics of NGC 7662 based on long-slit echelle 
spectroscopic observations and on Hubble Space Telescope 
(HST) archival narrow-band images. They inferred that PN 
7662 consists of a central cavity surrounded by two concen- 
tric shells, viz. the shocked dense shell and the outer slow 
AGB wind envelope, and estimated the temperatures of the 
two shells as ~ 1.4 x 10* K and ~ 1.1 x 10* K, respectively. 
By fitting the radial density profiles of the dense shell zone 
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Figure 5. Observed X-ray surface brightness profile (circles with 
error bars and our model best-fit curve (dotted curve) for the 
PN NGC 40 are shown here. The data was taken from Chandra 
X-Ray Observatory (ObsID 4480) in Montez et al. (2005). The 
radial profile was constructed by (i) regarding the X-ray centroid 
as the centre of NGC 40, (ii) summing X-ray photon counts in 
each of a specified set of concentric circular annuli, (iii) taking 
the error in each annulus for Poisson statistics (square root of the 
number of photon counts) , and (iv) dividing the counts and errors 
by the area of each annulus bin. The bin width for each annulus 
is chosen as ^ 3". A uniform background X-ray fiux ~ 3.4 X 10^^ 
ergs arcsec"^ is estimated in the circular region of radius 
> 25". 



Figure 6. Contours of the smoothed Chandra X-ray image over- 
laid on the optical BVR composite gray-scale image of NGC 
40 (Montez et al. 2005) with model-fitted dash-dotted R^s cir- 
cle (reverse shock) and dotted Rc circle (contact discontinuity; 
Rc > Rrs)- The X-ray and optical images are obtained by Chan- 
dra and WIYN observatory, respectively. The WIYN Observatory 
is a joint facility of the University of Wisconsin-Madison, Indiana 
University, Yale University, and the National Optical Astronomy 
Observatories. 



(see figure 13 of lLZ I ) and the slow AGB wi nd en velope with 
the self-si milar dynamic evolution model of ILZ I for isother- 
mal voids, ILZ I estimated relevant physical parameters for 
PN NGC 7662 which appear consist ent w ith optical obser- 
vations (e.g., Guerrero et al. 2004). ILZ I have also inferred 
properties of the inner fast wind by imposing condition (|20p 
for the pressure balance across the contact discontinuity. 
Shocked downstream inner wind velocity Vw,d = 25.9 km s~^ 
was inferred by the outer dense shell expansion. The inner- 
most fast wind mass loss rate M/u, « 2 x 1O~*M0 yr~^, the 
inner fast wind speed Vw,u = 1500 km s~^, the hot bubble 
(i.e. shocked downstream inner wind zone) plasma tempera- 
ture ~ 6.4 X 10^ K and the reverse shock velocity Vrs = — 10 
km s~^ represent a plausible set of fitting parameters that 
can reproduce a PN with more or less the same outer zone of 
NGC 7662. At the current stage of evolution, we have esti- 
mated the two radii Rc — 5.1 x 10^® cm and Rrs ~ 6.5 x 10^^ 
cm, i.e. the radius ratio Rrs/ Rc = 0.127 is fairly small. Equa- 
tion (g} gives the jS value as ~ 4.43 x 1O""M0 yr"V(km 
s~^). The total X-ray emissivity is A « 5.1 x 10"^'' ergs cm'^ 
s~^ for a plasma of temperature 6.4 x 10® K (e.g., Sarazin 
1986). For m^j = rap, we then have an X-ray parameter 
X = 1.42 X 10**® ergs s~^ cm. This would give a total X-ray 
luminosity of NGC 7662 as Lx = 2.4 x 10^^ ergs s~^ In 
short, we predict that an X-ray surface brightness like the 
upper-left plot of Fig. [1] with no distinct reverse shock sur- 
face and a total X-ray luminosity similar to that of NGC 40 
is anticipated for NGC 7662, which should be detectable by 
the Chandra X-Ray Observatory in space. 

However without available X-ray observations of NGC 
7662, physical conditions of its inner wind zone confined 
within the contact discontinuity radius Rc of the PN NGC 
7662 cannot be uniquely determined at this stage of model 



analysis. In other words, in terms of fitting available obser- 
vations of NGC 7662 in other bands (e.g. optical, infrared 
bands and so forth) , we can predictively explore X-ray emis- 
sions from the inner shocked downstream wind plasma of 
NGC 7662. In the above procedure, the innermost fast wind 
mass loss rate is a key yet adjustable parameter and may 
vary within a sensible range. We now show the correspond- 
ing ranges of our m odel results for NGC 7662 in Figs.[7]and 
[8] In fact, the lLZ I ISSV model solution with a parameter set 
{vrs = -10kms"\T„,d = 6.4x10'^ K, M/„ = 2x10"* Mq 
yr~^, v.w,-u, = 1500 km s~^} is only one of the plausible inner 
shocked wind solutions that meet pressure balance condi- 
tion (|20p across the contact discontinuity. By systematically 
varying the inner downstream shocked wind temperature 
Tm,tj and the reverse shock speed Vrs, we readily derive dif- 
ferent innermost fast wind velocities by equation ^ 
and the mass loss rates M/™ of the central star by wind 
pressure balance condition (|20|) across the contact disconti- 
nuity as shown in Fig. [T] To be specific, we pick a solution 
example marked by the plus symbol -I- in Figure [7] 

Figure [8] presents the upper limit of X-ray lu- 
minosity ~ M/u,u^ „, the predicted X-ray luminosity 
Lx by equation (|22|l and the ratio between these 
two quantities. The PSPC ROSAT count rate limit 
is ~ 0.03 s"\ If NGC 7662 is not detected in the 
ROSAT AU Sky Survey (RASS), then webPIMMS 
(http:/ /heasarc. gsfc.nasa.gov/Tools/w3pimms. html ) puts 
the upper limit for the X-ray luminosity of ~ 2.2 x 10^^ erg 
s~^ (a typical hydrogen column number density nn ~ lO'^^ 
cm"'^ and a hot bubble temperature Tra^d = 2 x 10® K 
are assumed for NGC 7662 at a distance of D = 1 kpc; 
e.g. Guerrero et al. 2004). The heavy solid curve in Panel 
B of Figure [8] marks this upper limit and suggests a 
possible regime for both Tru,u and Vrs at the upper right 
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Figure 7. Piecewise isothermal model solution properties of in- 
ner wind zone are explored for the PN NGC 7662. Panel A shows 
a chosen set of contour curves for logj^Q My„ in the unit of Mq 
yr~^ in a Vra versus T„ ^ presentation. Panel B shows a chosen set 
of contour curves for the inner fast wind speed Via^u in the unit of 
1000 km in a Vra versus T^, ^ presentation. The ^ axes of 
both panels are in the logarithmic scale. As discussed in the main 
text, the example model solution with parameters Vra = — 10 km 
s"\ T„,d = 6.4 X 10^ K, Mf^ = 2 X IQ-^Mq yr"!, = 1500 

km s~^ is marked here by the plus symbol +. 

side of the curve. In the regime where T^^u > 10^ K and 
Vra > —10 km s~^, the ratio Lx / {Mf^v%^^) is less than 
~ 10"^ = 1%. This estimate is also valid for PN NGC 40. 
Comparing Figures [7] and [8l we can also figure out that the 
mass loss rate A//„ < 10^^'^ M© yr~^ and the fast wind 
velocity Vw,u > 10'' km s~^. These estimated results appear 
all consistent with previous observations and theoretical 
calculations. 
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Figure 8. Energy solutions of inner wind zone properties for 
the PN NGC 7662. Following Figure^ here we show; (Panel A) 
a chosen set of contour curves for log^Q (Mj^„,vJ ^) in the unit of 
erg s~^ in a Vra versus presentation; (Panel B) a chosen set 

of contour curves for log^Q Lx in the unit of erg s"'^ in a Vra ver- 
sus Tu, li presentation; (Panel C) a chosen set of contour curves 
for logiglLx / {M f^v^ ^)] in a Vra versus T„ ,j presentation. For 
computing Lx , a minimum Rrs = 10 AU is adopted. The heavy 
solid curve in Panel B marks the upper limit given by ROSAT 
observation. The bottom left region to the heavy solid curve gives 
X— ray luminosity higher than the upper limit, which is forbid- 
den from the energetic consideration. The special model solution 
= -10 km s-i, T^^d = 6.4 X 10"^ K, Mf^ = 2 x W'^Mq 
yr~^, Vru^u = 1500 km s~^ is here marked by the plus symbol -|-. 



5 SUMMARY AND DISCUSSION 

We invo ke th e recently constructed ISSV hydrodynamic 
model of ILZ I to explore the dynamic shock interaction of 
the inner isothermal fast wind with the outer AGB slow 
dense wind and to calculate diffuse X-ray emissions from 
grossly spherical PNe for available observational compar- 
isons. Our piecewise isothermal wind dynamic interaction 
model involves a forward shock into the outer AGB slow 
wind envelope, a reverse shock into the inner fast tenu- 
ous wind and a contact discontinuity in between these two 
shocks. Different from previ ous dynamic model like Cheva- 
lier & Imamura (1983), |lz] model indicates that the X-ray 
emission peaks near the reverse shock, instead of the contact 
discontinuity surface between the hot bubble and the PN 
dense shell. Based on Chevalier & Imamura model, Akashi 
et al. (2006) indicates that diffuse X-ray emissions from hot 
bubbles should be in the form of narrow bright rings. Our 
piecewise isothermal model calculations show that the mor- 
phology of diffuse X-ray emission is determined by the ra- 



dius ratio Rra/Rc such that either ring- like or central bright 
X-ray morphology can appear. According to our dynamic 
model analysis, a relatively small Rra/Rc ratio is suggested 
for PN BD -(-30°3639, where a centrally bright diffuse X-ray 
emission morphology persists (e.g. Leahy et al. 2000; Kast- 
ner et al. 2000). In contrast, an adiabatic evolution model of 
PN cannot give a morphology for central bright diffuse X- 
ray emissions. We would expect various morphological cases 
for a more systematic survey of PNe. 

As current X-ray detectors in space have been powerful 
enough to observe extended X-ray morphologies of PNe, we 
offer quantitative line-of-sight X-ray brightness radial pro- 
files for candidate PN sources. In our model framework, the 
X-ray parameter X (defined in the main text), the reverse 
shock radius Rra and the contact discontinuity radius Rc to- 
gether completely determine both the X-ray luminosity and 
the surface brightness profile. In perspective, X-ray obser- 
vations together with optical and infrared images of grossly 
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spherical PNe are expected to give specific measurements or 
estimates for these three key parameters. 

The PN NGC 40 radiates intensely in UV bands, giv- 
ing a valuable observational platform to test and calibrate 
our X-ray model results for diagnostics. The results of our 
piecewise isothermal dynamic wind interaction shock model 
are fitted well with observations except that the central X- 
ray brightness is somewhat higher than that actually ob- 
served as shown in Fig. [5] NGC 40 appears to have a radius 
ratio Rrs/Rc close to unity with a relatively low temper- 
ature of ~ 10® K in the inner shocked downstream wind 
plasma. These two facts result in a fairly low X-ray lumi- 
nosity of ~ 4.0 X 10^^ ergs s"^ from NGC 40 (see Kastner 
et al. 2008 for this update). UV observations estimate a fast 
wind speed range of ~ 1800 - 2370 km s"^ (e.g., Cerruti- 
Sola & Perinotto 1985; Bianchi 1992). This is about twice 
our estimated fast wind speed ~ 1100 km s~^. We also note 
that the dipolar structure of NGC 40 significantly afi^ects the 
morphology of diffused X-ray emission. The actual hot bub- 
ble of NGC 40 is not very spherical, but in a dumbbell-like 
shape. Most X-ray photons are emitted from east and west 
regions, while only a few are from north, south or central 
regions. Observationally, spectral line analysis to the opti- 
cal bright rims show that the dense shell gas expands at 
~ 40 — 50 km s~^ in the north-south direction and at ~ 30 
km s~^ in the east-west direction (e.g. Sabbadin et al. 1999). 
This indicates that the earlier AGB slow wind is denser in 
the east-west direction so that the collision between the in- 
ner fast wind and the slow AGB wind is more violent in the 
east-west region. This also explains why the dense shell is 
brighter in the east and west regions. By both quantitative 
calculation on spherical smoothing and qualitative analy- 
sis on dipolar structure, we conclude that wind-wind shock 
dynamic interaction appears natural for the hot bubble for- 
mation and diffused X-ray emissions. 

Different from the PN NGC 40, the PN NGC 7662 is 
here inferred to have an inward reverse shock velocity of 
~ -10 km s"^ By expression (gU, NGC 7662 is predicted 
to have an increase of X-ray luminosity with increasing time 
t at an estimated rate of ~ 0.9 x 10^** erg s~^ per year. For 
example, its X-ray luminosity will reach ~ 2.1 x 10^*^ ergs 
s~^ 20 years later and ~ 2.5 x lO^'' ergs s~^ 50 years later. 
Energy conservation shows a minimum Rrs exists so that 
Lx must be limited by a maximum ~ M/^u^ ^. In fact, 
studies on NGC 40 and NGC 7662 both lead to the conclu- 
sion that Lx < 10~^A//u,i;2),u- The sensible explanation is 
that cooling effect becomes very significant even when Lx 
is still far lower than this upper limit. 

At the present epoch, it appears justifiable to not in- 
clude cooling effects in our piecewise isothermal dynamic 
model calculations. Akashi et al. (2006) considered the ra- 
diative cooling semi-quantitatively for computing the X-ray 
luminosity based on the self-similar colliding adiabatic wind 
model of Chevalier & Imamura (1983) and showed that the 
radiative cooling was unimportant during most time of a 
PN dynamic evolution. Nevertheless, as radiative cooling be- 
comes important for the hot bubble evolution, e.g. ~ 10'^ yrs 
later for the case of Vrs = —3 km s^^ in Figure[3l it would be 
inaccurate to compute X-ray emissions either in our model 
or in that of Akashi et al. (2006). The X-ray hot bubble 
will cool down much faster, the contact discontinuity bal- 



ance condition would be broken and the outer dense shell 
dynamic behaviour may be also interrupted. 

Since X-ray luminosity only takes up < 1% of M/^n^^,, 
the energy conservation consideration in Section 3 implies 
that a significant amount of energy should have been con- 
verted to other forms during the entire life of PNe. Oth- 
erwise, the temperature of shocked downstream inner wind 
should be ~ 10 times higher than the observed tempera- 
ture inferred from diffuse X-ray emissions. Ifowever, rela- 
tively low temperatures of X-ray emitting gas are gener- 
ally inferred in PNe at different evolution stages. We note 
that the energy conservation calculation in Section 3 is 
fairly general, where we focus on the mass conservation 
Mf^ = 'i-Kr^ pw,dVw,d and choose typical observed param- 
eter values without specifically invoking either the isother- 
mal or adiabatic assumptions. Therefore we suspect that 
there must be certain long lasting mechanisms to explain 
the low temperature of X-ray emitting gas, and such mech- 
anisms are not quantitatively taken into account in either 
[lz] piecewise isothermal model or the adiabatic model of 
Chevalier & Imamura (1983). Moreover, such mechanisms 
influence more than radiative cooling does so that shocked 
hot bubbles do not evolve adiabatically, as implied by the 
observation of central bright X-ray emission observation in 
PN BD +30°3639. Thermal conduction and material diffu- 
sion might be non-negligible if the magnetic fleld near the 
expanding contact discontinuity surface becomes sufficiently 
weak. Another possibility is that the hot bubble largely be- 
tween the reverse shock and the contact discontinuity be- 
comes turbulent to partially account for the missing energy. 
The thermal conduction between hot bubble and cool dense 
shell could be considerable (e.g. Kastner et al. 2000; Steffen 
et al. 2008). The rapid variance of inner fast wind mass loss 
rate of central stars might also be essential to account for the 
imbalance between fast wind kinetic energy input and low 
temperature inferred by diffuse X-ray emissions (e.g. Soker 
& Kastner 2003; Akashi et al. 2007). In Sections 2 and 3, we 
did not consider such mechanisms; instead, we choose the 
shocked fast wind temperature Tw,d (or X) as an adjustable 
parameter for fitting X-ray observations. 

In some PNe with detected X-ray emissions, very dis- 
tinct bipolar morphology appears, such as HENIZE 3-1475 
(e.g. Sahai et al. 2003), NGC 7026 (e.g. Gruendl et al. 2006), 
MENZEL 3 (e.g. Kastner et a. 2003), and so forth. CoUi- 
mated outfiows and jets in magnetized gas media might be 
expected to account for these types of diffuse X-ray emis- 
sions. Akashi et al. (2008) performed two-dimensional nu- 
merical simulations of jets expanding into the slow wind 
of AGB stars. They proposed that this jet-wind or outflow- 
wind interaction model might explain such bipolar morphol- 
ogy of diffuse X-ray emissions. 
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